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Abstract 
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1 Introduction 


Pulsar timing is a promising method for the detection of stochastic back¬ 
ground of gravitational waves (GWs) (see, e.g., Hobbs et al. [lj, Wang et 
al. p]). As shown by Hellings and Downs [3j, the shift in the pulses arrival 
time can be related to the spectrum of this background. Pulsars also exhibit, 
besides very regular pulses, important properties such a the polarization of 
their signals. This polarization has a complicated structure (see, e.g., Craig 
m, but it can be affected by GWs and it may be possible therefore to ex¬ 
tract some information about their presence in the galaxy. As shown in the 
present paper, this could be achieved by adapting and complementing the 
statistical analysis proposed by Hellings and Downs [3j. 

The aim of this paper is to analyze the interaction of electromagnetic 
waves with a plane fronted GW, and to deduce the corresponding formulas 
for the Stokes parameters and the rotation angle of the polarization. A 
similar problem was studied from a formal point of view in a previous paper 
[5], with the GW described by the Ehlers-Kundt metric [6], which is an exact 
solution of the Einstein equations. It was shown that there is a rotation of 
the polarization angle in phase with the GW. In the present paper, we apply 
a similar analysis to a GW described in the standard form with its two 
independent polarization states. Section 2 is devoted to the mathematical 
formalism in the short wave-length approximation, and the electromagnetic 
coherency matrix is obtained in Section 3. A possible application to the 
detection of stochastic GWs background is discussed in Section 4. 


2 The electromagnetic field 


The metric of a plane gravitational wave in the weak field limit is 
ds 2 = —2 du dv + (1 + f)dx 2 + (1 — f)dy 2 + 2 g dx dy. 

The relation with Minkowski coordinates t and z is 

1 , 1 , 

u = -j=(t- z ), v=-=(t + z), 


( 2 . 1 ) 


and the two degrees of polarization of the GW are given by the potentials 
f{u) and g(u), which are functions of u only. In the following, quadratic and 
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higher terms in / and g are neglected, coordinates are taken in the order 
(u, v, x, y ), and we set c— 1. 

In the short wave-length approximation (see, e.g., Misner, Thorne and 
Wheeler [7J or Ref. 0). the electromagnetic potential four-vector A a is 
written in the form 

A a = a a e tS , (2.2) 


in terms of an auxiliary four-vector a a and the eikonal function S which 
satisfies the equation 

g^S^S, u = 0. (2.3) 


Defining the null four-vector 


AT 


S 

? 


cn 


the Maxwell equations can be expanded in descending powers of S. The first 
term in the expansion implies 


a a K a = 0 (2.4) 

and the electromagnetic field turns out to be 

F afi = i(K a ap - Kpa a )e lS . (2.5) 

The second term in the expansion implies 

K l3 (a a ,p ~ ap, a ) + ( a a K 0 - a /3 A' Q ) ;/3 = 0. (2.6) 


Choosing a gauge A v = 0 as if Ref. [5], a particular solution of the above 
equations turns out to be 




1 

kn, - 


a v = 0 


kx^X kyCiy 

1 


d x = (l + -/(«)) d x + -g(u) d y 
a y = kg{u) d x + (l- ^/(u)) d y 


where d x and d y are constants. 


As for the null-vector K a , its components are 

AT •A('u)] , k v , k x , kyj , 


(2.7) 


( 2 . 8 ) 
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where 


J-’(w) = f(u) cos 20 + g(u) sin 20, (2.9) 

k a = (uj, k) is the wave-vector in Minkowski coordinates, and the angle 0 is 
defined by k x = k± cos 0 and k y = k± sin 0, with k± = k x + k y . In particular 

, oj + k z 

u ~ V2 

j u ~ k z 

v ~ V2 ' 

(In the limit of flat space-time, the above solution describes a plane elec¬ 
tromagnetic wave with the electric field 

E = z(caa — a 2 k)e*^, (2.10) 


where S = k Q x a = —c ot + k • r and oja z = k ■ a (see the Appendix of Ref. 

0)0 

We can now choose a time-like unitary vector 


t“ = -^[l,1,0,0 


( 2 . 11 ) 


which is also a geodesic, that is, it is tangent to the world-line of an inertial 
observer fixed at a point with constant coordinates (x,y,z). 

The observed frequency of the electromagnetic wave is defined covariantly 


as 


fi(w) = ~K a t a = —={K U + k v ) = 
v 2 


o;[l - -(1 + cos0)J 7 (m)], 


( 2 . 12 ) 


where 6 is the angle between k and the direction of propagation of the GW, 
and of course to = —k a t a , which is the frequency observed in the absence 
of the GW. From this last formula, it follows that the Doppler shift for a 
pulsar located at a distance l from Earth and at an angle 9 with respect to 
the direction of propagation of the GW is 


[0(f) — 0(f — l — l cos 9)\/uj, 

which is the formula deduced by Estabrook and Wahlquist [8] for the pulse 
frequency of a pulsar (for / type polarization only). 
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3 Coherency matrix 

Given the time-like vector t a , the electric held can be defined as E a = F a/3 tp, 
where F Q/3 is the electromagnetic tensor. In the short wave-length approxi¬ 
mation [5j 

E a = i[Qa a + (apt l3 )K Q \e tS . (3.13) 

Two space-like unitary vector e«' ) (i = 1,2), orthogonal to each other and 
to both t a and K a , can be dehned. Accordingly 

A Q e (i)Q = ifia a e (i) “e iS . (3.14) 


A possible choice satisfying all the orthogonality conditions is 

1 


e (1)Q = 


k_ L\/l — T 


0, 0, kyj 


and 


el 2 ) — — -tt—— r— 1 — 1,1,0,0 


K u + k v 
1 (. K u - k v ) 


k±s /1 — T K u + k v 


0 , 0 , kxi ky j . 


Accordingly, we have 


a„e (1) “ = 


k±y/l — T 


( &xky + Ojy kx') 


( 2 )a _ a uk v 


a„e' = 


Wi — F 


and therefore 


£ (1) = E n e {l)a = i 


Q 


k_ l\/1 — F 


( Clyk X dxky^jC 


iS 


E {2) = E n e^ 2)a = i- 


n 


-.a u k v e 


iS 


k ± Vl — F 

These formulas can be written in the compact form 


(EW\ ,^e iS r 1 

(s«) - ! TT [ + 2 6 


0 1 

-1 0 


fox \ ( ^.1 


fcyJ \^yJ 


(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 
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where 


(3.22) 


Q{u) = g{u) cos 20 — f{u) sin 20, 
or even more compactly as 


(E^ 

\ew 


= i-jV-S) 

LU 



G 





(3.23) 


where E W are the values of the electric field in the absence of the GW. Thus 
the contribution of the GW is made implicit through the functions E{u) and 
G{u), and it follows from this last formula that the GW produces a rotation 
of the polarization vector by an angle —(1/2 )Q(u). 

The coherency matrix is 

8 - ( ^ •% |:f) . (3.24) 


where S) are the Stokes parameters (in the notation of Born and Wolf [9]). 
It then follows with some algebraic manipulations (such as in Ref. [5]) that 

S= S § - 6 (t 2 |). ( 3 ' 25 > 


where § is the coherency matrix and S) the Stokes parameters in the absence 
of the GW. 

More explicitly 


S 0 = So[l - (1 Tcos^J 7 ] (3.26) 

Si+iS 2 = (Si + iS 2 )[l - (1 + cos6)E - iQ\ 

S 3 = S 3 [l - (1 + cos#)^ 7 ]. 

It can be seen from these last formulas that if the light is not initially 
polarized, namely < S) >= 0 for i = 1, 2, 3, it will not be polarized by the 
GW. However, if the light is originally polarized, the angle of the polarization 
ellipse oscillates at the passage of the GW. This can be seen by setting 

S 1 + iS 2 = S 0 cos2y e 2 ^ 

S 3 = S 0 sin 2x, 
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in terms of the Poincare angles (see, e.g. 0): 0 specifies the orientation of 
the polarization ellipse and x its ellipticity. It thus follows that the ellipticity 
of the polarized light is not altered by a GW, but the polarization angle 0 
changes as 

0^0-^(u), (3.27) 

as already implied by (13.2311 . 

The intensity of the wave changes as So and the variation of the frequency 
is given by (12.12(1 . 

4 Correlations 

It is possible to adapt the analysis of Hcllings and Downs p] to the polar¬ 
ization shift as given by the above formulas. The change in the polarization 
angle A0 observed on pulsar number i is 

A 0i = -ong(t) + Pif(t ) + 7100, (4.28) 

where f(t) and g(t) are the GW signals common to all pulsars, 

OLi = ^ COS 20j, 

A = ~ sin 20j, 

are angles factors for the i th pulsar, and n00 represents all noise sources 
unique to the pulsar, including terms depending on its position. 

The corresponding analysis is just as in Hellings and Downs p). The 
correlation < A0*A0j > for two pulsars is 

Cij = aiaj < f 2 > +Pifij < g 2 >, (4.29) 

in terms of the correlations of the potential functions, the basic assumptions 
being that the processes n,;, rij, together with f(t) and g(t), are all uncorre¬ 
lated among themselves. 

For an isotropic background of GWs, one has to calculate the angle av¬ 
erages for pairs of pulsars 

If If 

an = — / cqctq dSl = — 

4vr ./ 4 t t J 
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/3i/3j dSl. 


(4.30) 






The result is 


11, /I — cos7i,-\ 

a H = g + 4 I 1 - cos 7«j )ln ^- 2 - ) + 

-(1 + cos 77 ) In ^---J, (4.31) 

where 7 \j is the angle between the two pulsars (see the Appendix for an 
outline of the derivation). 

The rest of the formulation is mutatis mutandis just as in Ref. [3] and 
will not be repeated here. 

5 Conclusion 

The influence of a GW on the polarization of light is made explicit by the 
formulas obtained above. The basic effect of a GW is to produce an oscillation 
of the electromagnetic wave polarization. This effect could be applied to the 
detection of GWs. An example of such an application was outlined as a 
complement to the pulsar timing method. 


Appendix 

Consider first gravitational waves propagating along a particular direction, 
say the z axis. An average of the product a±a 2 over the azimuthal angle 0 is 

1 f 2rc 1 

— / dcf) cos( 2 ( 0 i - <j>)) cos(2(0 2 - </>)) = - cos( 2 ( 0 i - 0 2 )), 

07T J q o 

where are the azimuthal angles of the pulsars. The same result applies 
to the averaged product /3 i/ 3 2 - The above expression can be written in a 
coordinate independent form as 

[cos 7 - (a! ■ g)(a 2 ■ g )] 2 _ 1 

4[1 - (ai • g) 2 ][1 - (a 2 • g) 2 ] 8 ’ 

where a,; are the unit vectors in the directions of the pulsars and g is the unit 
vector in the direction of the considered GWs. 

If 7 is the angle between the two pulsars, we can now choose without 
losing generality: 


a : = (cos( 7 / 2 ),sin( 7 / 2 ), 0 ) 





a 2 = (cos( 7 / 2 ),-sin( 7 / 2 ), 0 ), 

and it then follows that an integration over all the directions g yields 


Ot 12 — 


167T 


dVt e 


(cos 7 - g^) 2 

L(1 ~ ~ 9 %) 


1 
2J ’ 


where 


gi = sin# cos (0 — 7 / 2 ) 
g 2 = sin 9 cos+ 7 / 2 ) 

Using as an intermediary step the general formula 


f2vr 


df 


cos 7 — ( 1 / 2 ) sin 2 9 (cos( 20 ) + cos 7 ) 


1 2 


[1 — sin 2 9 cos 2 (0 + 7 / 2 )] [1 — sin 2 9 cos 2 (0 — 7 / 2 )] 
cos#! (1 + cos 2 9) 


2vr 


1-2 


sin 4 9 + 4 cos 2 0 / sin 2 7 


one obtains Eq. f!4.3ip . 
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